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Abstract The Heisenberg-Euler Lagrangian (HEL) is not only a topic of 
fundamental interest, but also has a rich variety of diverse applications 
in astrophysics, nonlinear optics and elementary particle physics etc. 
We discuss the series representation of this Lagrangian and a few of its 
applications in this study. [In an appendix, we discuss issues related to 
the renormalization — and the renormalization-group invariance — of 
the HEL and its two-loop generalization.] 
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1 INTRODUCTION 



The effective one- loop Lagrangian density [1-11] of quantum electrodynamics (QED) describes 
the nonlinear interaction of electromagnetic fields due to a single closed electron loop. [It is also 
referred to as the QED effective action, or the effective Lagrangian.] One of the most compact 
and elegant ways to treat the symmetry properties of the vacuum is by the method of the 
effective Lagrangian 12 . 

This one-loop Lagrangian, often called the Heisenberg-Euler Lagrangian (HEL), has been used 
to describe a variety of electromagnetic processes. The real part of the Lagrangian can be 
used to delineate such dispersive phenomena as photon propagation in a magnetic field, second 
harmonic generation, photon splitting in a magnetic field, and light scattering in a vacuum 
[1,4-9,11,13-17]. The development of high intensity lasers, with the consequent availability of 
powerful, coherent light sources, can render possible the observation of delicate nonlinear effects 
and also the nonlinearity of Maxwell's equations. Strong magnetic fields around pulsars and 
around magnetars, and other astrophysical and laboratory situations might also warrant accurate 
calculations based on the Lagrangian for astrophysical and other applications. The imaginary 
part of the Lagrangian has been applied to absorption processes such as electron-positron pair 
creation gllSllHlin] • 

While the integral representation of this Lagrangian as a function of the electric and magnetic 
field is well known, further analytical representations and numerical procedures for arbitrary field 
strengths have recently been refined. In the special cases when either the magnetic or electric field 
vanishes, the Heisenberg-Euler Lagrangian may be expressed in terms of elementary functions 
and integrals of the natural logarithm of the generalized gamma function with real or complex 
arguments |18ll9j . Elementary function series expansions for these gamma function integrals and 
precise numerical values for the Lagrangian using the expressions of Dittrich et al. were provided 
by Valluri et al. |[lSj. For the more general case when both an electric and magnetic field are 
present, a numerically useful, analytical series expression for the real part of the Lagrangian was 
derived and this latter series expression involves elementary functions and the sine, cosine, and 
exponential integrals, all of which are easily calculated |2U| . 

Motivations for studying non-linear generalizations of Maxwell's equations in the vacuum are now 
quite different. Strong interest in one-loop corrections to the classical Lagrangian in Abelian 
(as well as non- Abelian) gauge theories has resurged in order to learn about the structure of the 
vacuum when it is probed by an external electromagnetic field. The problem of the existence 
of a stable electron has been and continues to be an interesting issue, not only relevant to pure 
electrodynamics but also to an unknown fundamental theory of matter and its interactions (for 
an interesting view on this issue, see |12j). 

In this paper we discuss some analytic calculations relevant for the applications of the HEL. The 
paper is organized as follows. Section 2 provides the analytical definitions and preliminaries. 
Section 3 presents the analytical results indicating the existence of the higher harmonics for 
the real part of the HEL, and the discussion of these results will conclude the main body of 
the paper. Some aspects of the renormalization and the renormalization-group invariance of the 
action are briefly discussed in an Appendix. The final section summarizes the conclusions. 
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2 REPRESENTATION OF THE QED EFFECTIVE 
ACTION BY SPECIAL FUNCTIONS 



The renormalized (see also the Appendix) Heisenberg-Euler Lagrangian (HEL) IS.C is expressed 
as a one-dimensional proper-time integral [1-4]: 

■ , w ^ . , ^ a' - 1 1 

ao coth(eas) cot (eos) — (1) 

3 (es) J 

and is a quantum correction to the Maxwellian Lagrangian 

La = Lo = Lm = \{E^- B^) ■ (2) 

A>C is a one-dimensional proper time integral with proper-time parameter s and can be evalu- 
ated by numerical quadrature. A convergent series expansion with (computational and 
conceptual advantages) can be derived in terms of special functions. There is a unified series 
expansion encompassing both the real and the imaginary part [^. We need techniques for a 
reliable numerical evaluation which may have a variety of applications. While the series 
expansion might be regarded as a complete solution from a theoretical point of view, the expan- 
sion alone does not solve the physics associated with the Lagrangian, and it does not provide eo 
ipso a general, numerically efficient algorithm for its evaluation. The integral defining the HEL, 
which can be calculated numerically by quadrature, has a slow convergence due to oscillations 
of the integrand. 

The QED effective Lagrangian can be expressed as a function of the Lorentz invariants J- and 
Q, given by 

T = \f,^F^^ = \^^"- = ^ - b') , (3) 



A£=-At lim / ^^-{m^-^e)s 



G = \f^u {*FY'' = -E-B = ±ab, (4) 

where E and B are the electric and magnetic field strengths, F^i, is the field-strength tensor and 
{*Fy denotes the dual field-strength {*FY'' = {l/2)e'"'P'' Fp„. The quantities a and b denote 
secular invariants, 

a = \/Vj'^ + G^ + :F, (5) 

h = ^ yjj^ + g^-j^. (6) 

Secular invariants emerge naturally as eigenvalues of the field-strength tensor fields. We also 
note that in the case ^ < 0, it is possible to choose a Lorentz frame in which E and B are 
parallel. 

a = \B\ and h = \E\ if B is (anti-)parallel to E , (7) 
a and b are positive definite and this gives the condition: 

ab = \E ■ B\ > {) for any Lorentz frame and Q ^ Q . (8) 
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Also, the following notation is sometimes used: 

T = -s, g = -p, 



(9) 



= -^= -If^.f^" = ^ - B") = ^ - • (10) 

The correction AC to the Maxwellian Lagrangian Lq can be written in terms of the secular 
invariants a and b [11121) 

AC = ReAC + ilmAC . (11) 

the real part can be expressed as 

2 ^ 

Re AC = -^abY^ [an + dn] , (12) 

n=l 

coth(n7r6/a) ( „. f mTm'^\ / nvrm^X . / nvrm^X . / rmrm 



n 



{ ^ f mrm \ ( mrm \ ( mrm \ ( rmrm \ 1 , , 

Ci cos + si sin , 13 

I \ ea J \ ea J \ ea J \ ea J \ 



— coth(n7ra/6) I f n-Km?\ „ ( n'Km?'\ 



2n 



exp 



n-Krr? \ „. / rmrrr? 



eh 



while the imaginary part is given as 



n=l ^ ^ 



(15) 



We note that the Cosine and Sine integrals have an "asymmetric" form since the generally 
accepted definitions for these integrals are "asymmetric" and are shown below. More details are 
given in |20j . 



oo 



Ci(z) = -/ dt^^ z>0, (16) 



oo 



si(z) = -/ dt^ = Si(z)-|, (17) 

Si(.) = Tdt'-^. (18) 
Jo t 

The imaginary part ImAiZ is generated by a modification of the integration contour in the 
exponential integral entering into the definition of d„ (normally, the exponential integral is 
defined via a principal- value prescription). The relevant exponential integral is 

— ; — in the definition of dn and reads (19) 
eb ) 

/■oo -t 

Ei(u) = -(P.V.) / dt forueR. (20) 

.1 —oi t 
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Under an appropriate deformation of the contour, a unified representation for both the real and 
the imaginary parts is obtained j2()j . 

2 °° 

A/: = hm V [6„ + c„] , (21) 

n=l 

coth(n7r6/a 

On — 



2n 




(22) 



coth (mra/b) 
"""" 2^^^ 

(23) 

The effective action has branch cuts along the positive and negative 6-axis as well as along 
the positive and negative imaginary axis. For more details, the reader is referred to |2U| . In 
contrast to the exponential integral Ei, the incomplete Gamma function is defined in the entire 
complex plane with a cut along the negative real axis. It is important to see this connection 
since the Barnes function |21j is closely connected to the Gamma function which has numerous 
applications. 

The main numerical difficulty is the slow overall convergence of the series expansion whose terms 
are of nonaltering sign. Fade approximants, a standard tool in many power series application 
are not capable of summing the series for the HEL |2nj . The terms of the convergent series 
representation are interpreted as being generated by a "partial-fraction decomposition" 
in close analogy to "partial-wave decomposition" in bound state calculations. It has been 
shown that the convergence of the HEL series can be accelerated by the same technique as in 
partial wave decomposition also called the "combined nonlinear-condensation transformation" 

(CNCT) [laEaEni- 

Another interesting representation is in terms of the Barnes function. The Barnes function G(z) 
is a generalization of the Euler gamma function and is also related to the Hurwitz Zeta function. 
It obeys the recursion relation [2j 



G(z + 1) = F(z)G(z) zee, G(l) = l. (24) 
We write a useful expression for logG(2; + 1) below: 

log G{z + 1) = log(2vr) - 2|! _ '±^1 

+ (25) 

k=2 

logG(z + l) = zlog(F(z)) + C'(-l)-C'(-l,^), (26) 
C{t,z) = j^C{t,z), (27) 
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where C{k) is the Riemann Zeta function, used in the representation of the HEL. Due to its 
close connection to the T{z) and the Zeta function, the Barnes function may have a variety of 
apphcations in Computer Algebra and Theoretical Physics as well as in other fields. 

The Mittag-Leffler theorem (see PUj for a comprehensive discussion) is a key element in the 
derivation of the series representation for the HEL (some useful formulas, originally derived 
without the explicit use of the Mittag-Leffler theorem, still based on a partial-fraction decom- 
position and equivalent to the results obtained by us using the Mittag-Leffler theorem, can be 
found in \2i[ p. 271]). The generalization of the Mittag-Leffler theorem for a function with non- 
simple poles and zeros is the Riemann- Roch theorem j25(l26j . For Lagrangians with poles and 
zeros of higher order, which can occur in string theories, we anticipate the application of the 
Riemann-Roch theorem. Such Lagrangians with poles and zeros of arbitrary order may form a 
vector space of meromorphic functions over a complex field. 



MAGNETO-OPTICAL EFFECT: SECOND 
AND HIGHER-HARMONIC GENERATION 



Maxwell's equations receive corrections from virtual excitations of the charged quantum 
fields(notably electrons and positrons). This leads to interesting effects PP: light-by-light scatter- 
ing, photon splitting, modification of the speed of light in the presence of strong electromagnetic 
fields, and - last, but not least - pair production. The dominant effect for electromagnetic fields 
that vary slowly with respect to the Compton wavelength (frequencies u! <C 2m(? /h) is described 
by the Heisenberg-Euler Lagrangian, which is known to all orders in the electromagnetic field. 
For the case of zero electric field, the HEL can be written as jl8lll9j 



A/:(/i) 



m 

32^ 



- + 2h + 2h' 
3 



In/i + - 4Li + 41nri (1 + h) 



(28) 



where /i = w = - 0-2487, H^r = 4.4 x lO^^ Qg^^gs and Ti{h) is the generahzed 

r function. We observe that the expression in the braces for the HEL contains terms with 
quadratic powers of h as well as terms like hlnh and /i^ In h. It is well known that the logarithm 
of a trigonometric Cosine function has even powers of the argument when expressed as a power 
series |161l27l05| . This is the situation encountered for AC{h) when the magnetic field H (also 
meant as B) is expressed as a cosine wave for a plane electromagnetic wave propagation jl6j . 
The cases of parallel and perpendicular polarizations for the weak and strong field cases and 
the evaluation of the higher harmonics have been treated in detail earlier in JH]- Equation (28) 
describes the case for arbitrary field strengths and generalizes the evaluation for the generation 
of the higher harmonics. 

Magneto-optical Effect We briefly discuss second and higher-harmonic generation in a static 
magnetic field which is an example of broken symmetry. By broken symmetry, we mean a 
violation of the principle of superposition/linearity. High intensity synchrotron radiation, storage 
rings and lasers are good experimental tools for the study of higher-harmonic generation. The 
degeneracy of nonlinearity is broken by a static field or nonplanar or nonmonochromatic waves 
and leads to the generation of higher harmonics. We define. 



1 



X 



1 



S + {S' 



(29) 
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in analogy to a and b of equations (5) and (6) above. Here Fc = Her- For 

X = 0, y 7^ 0, 1^1 = 0, we obtain 

TT ^ — ' 



(30) 



k=l 



where, 



ey Ei{ — + e y Ei [ — 

y J \y 



N=l 



(31) 



Again, it is to be noted that Ei in Equation (31) is defined as real for real (negative and/or 
positive) arguments. 

We use the formulas: C(4) = ("(6) = ^ etc. to derive the equations above and below, 
can be evaluated in the case of an external static magnetic field Bq along the x-axis. For an 
electromagnetic wave propagating along the z axis, we have for the case jT], 



Ex = Ez = 0, By = Bz = 0, \Ey\ 



\B^ 



p = E-B = 0, 



S = I [El-{Bx + B,) 



'2 (-^0 + 2-Bx-Bo) 



We then derive an expression for the partial derivative of the ReCeff with respect to S: 



1+^ 

vr 



4 5 16 52 256 ^3 



+ 



ASF^ 105 315 F6 

The electric displacement field D can then be obtained. 



+ 



d 

D ReCeff 



■E. 



We also observe that 



dP d 

^•^Re£.,, = 0. 
We consider electromagnetic waves with Ey and B.^ of the form: 



We then find that. 



Dy = Ey 



1 + 



Ey 

Bx 



4a 
457r 



A cos cj} = A cosi^wt — kz) , 
A cos (j) = A cos{wt — kz) . 



Bq +2Bx Bp 
2F2 



+ 7 



3 f Bq + 2 Bx Bq^ ^ 



F2 



(32) 
(33) 

(34) 

(35) 
(36) 

(37) 
(38) 

(39) 



where A is the amplitude of the wave. We find that the parallel mode propagates alone. The 
EyBx terms in equation (39) indicate the presence of the second harmonic, as can be observed 
from the fact that: 



2EyBx = 2A'^ cos^ ((> = A'^{1 + cos 2(f>) 



(40) 
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The term proportional to cos(2(/)) indicates the presence of the second harmonic. Similarly, EyB"^ 
indicates the presence of the third harmonic and EyB^ indicates the fourth harmonic etc. In a 
similar way, the evaluation of H can be done: 

^ = -||-^Re£e//(P = 0). (41) 

The above expression of H is also shown to exhibit the second and higher harmonics. The 
spatial anisotropy is imposed by the magnetic field. The SHG depends on propagation direction 
and polarization direction of the fundamental electromagnetic wave with respect to Bq. The 
maximum effect is when Bq is -Lk (the wave propagation direction). 



4 CONCLUSIONS 

We have investigated questions related to the representation of the quantum electro dynamic 
(QED) effective Lagrangian and its analytical expression. In Sect. 2, we briefly recalled our 
previous results given for special-function representations of the effective Lagrangian, and we 
briefly clarify the mathematical notation used in the special-function representations (12) and 
(21). The representation (21) unifies the real and imaginary parts. We also introduced the Barnes 
function that should be of use in a variety of applications that include the HEL. We very briefly 
mention that the key step in the derivation of our series representations is the Mittag-Leffler 
theorem j2Uj . 

In section 3, we have also given an expression from the effective Lagrangian that will facilitate 
the evaluation of the higher harmonics. We also briefly discuss the magneto-optical effect of 
the vacuum which might provide a signature of "QED's nonlinear light." In the Appendix, we 
discuss some aspects of the renormalization and the renorm-invariance of the action. Based 
on the results of the current paper, we expect to carry out detailed studies related to various 
projected and ongoing experiments and astrophysical phenomena [1,11-17] involving strong 
static-field conditions (or fields with frequencies that are small as compared to the electron 
Compton wavelength). The study of the HEL to include finite temperature effects |29) with the 
series representation that we have developed is an interesting problem that warrants further 
investigation. 
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A APPENDIX - RENORMALIZATION 

AND RENORMALIZATION-GROUP INVARIANCE 

The renormalization as well as the renormalization-group invariance of the QED effective La- 
grangian have already been discussed at length by various authors (see [2111 and references 
therein). 



8 



However, material on this topic is somewhat scattered in the literature. In the current Appendix, 
we would like to give a brief overview of the physical ideas that led to the formulation of the 
related RG equations, as well as the connection with the (3 functions of QED, complemented by 
an easy-to-understand presentation of the basic notion of the renormalization itself. 

It should be noted that the considerations in this section have only partial significance for the 
main topic of the current investigation, which is a treatment of higher-harmonic generation based 
on series representations of the Lagrangian. The current section merely provides illustrating 
remarks on the derivation of the effective Lagrangian, which is again the starting point of the 
main endeavour pursued. 

The "original" Lagrangian of QED from which we start to work out the S-matrix with all its 
radiative corrections in fact has to be identified with the "bare" Lagrangian that is expressed 
in terms of the "bare" physical quantities (bare charges, masses, fields). As we develop the 
perturbations series (in e^), we see that some terms (it does not matter if they are infinite or 
not) in fact modify the physical parameters that entered into the very Lagrangian from which 
our work started. This means that we have to renormalize the Lagrangian. 

Renormalization: The renormalization can be done by adding counterterms to the Lagrangian. 
These terms relate the bare parameters (bare charges, masses, fields) to the physical parameters 
of the theory, i.e. to the renormalized charges, masses, fields. The renormalizability then requires 
that the Renormalized Lagrangian = Bare Lagrangian + Counterterms = Bare Lagrangian but 
with charges, masses, fields expressed in terms of the renormalized parameters. For QED, we are 
in the lucky position that the theory remains finite in the infrared, i.e. that the physical charge 
of the electron remains finite as we move two electrons far apart. Therefore, we may renormalize 
QED on mass shell, i.e. renormalize QED in such a way that the renormalized charge is by 
definition the charge of an electron as seen by another electron in the limit of a large separation 
of the two electrons (this latter situation corresponds to the limit of a very soft exchanged photon 
— the infrared limit in which the wave vector k of the exchange photon tends to zero). 

To see how this works, look at p. 325 of The vacuum polarization modifies the Coulomb 
law (or Thompson scattering) in such a way that (eo is the bare charge) 

el/k^ - el/[k\l + uik-"))] = el/[k\l + u;{0) + Oik'))] . (42) 

Now at A; = 0, the physical charge of the electron is obtained, which means that = eg/(l + 
Lij(O)), or that the renormalization that relates the bare and the renormalized charge of the 
electron is obtained as = l/(l+a;(0)) where = Z^e^. This relates the physical renormalized 
charge with the bare charge of the electron. In order to construct a renormalized Lagrangian 
in which the "e" is indeed the renormalized, finite electron charge, we now have fulfill the 
condition that the vacuum polarization as derived from our renormalized Lagrangian results in 
a "modified" vacuum polarization correction to the Coulomb law that fulfills (^Renormalized (0) = 0. 
This is just Eq. (8-96e) of In general, a set of renormalization conditions imposed on the 
renormalized Lagrangian determines the physical interpretation of the parameters in terms of 
which the renormalized Lagrangian is written. The renormalization conditions (8-96) in |^, on 
which the considerations of |2Qj are based, ensure that the renormalized effective Lagrangian is 
written in terms of the finite, renormalized, physical electron charge. The "usual" representation 
of the QED Effective Lagrangian (QED Effective Action) in Eq. (1) fulfills these conditions. 

In other words, since we are used to defining the value of the coupling by Thompson scattering, 
i.e., by scattering a long- wavelength (uj 0) photon off a (static) electron, the renormalization 
scale is naturally set by the mass of the electron. This amounts to the so-called on — electron — 
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mass — shell renormalization condition e^(/i = m)/(47r) ~ 1/137.036. Our starting point Eq. (1) 
for the HEL is obtained exactly by implementing this renormalization condition, and the coupling 
e used throughout this work should be interpreted in this way. It is precisely the second term in 
square brackets in Eq. (1) that guarantees this renormalization condition. This term subtracts 
any contribution of AC to Thompson scattering; therefore, the latter is fully described only by 
the renormalized Maxwell term together with the static-eletron interaction term involving the 
on-shell renormalized coupling e - as it should. This renormalization procedure has thereby fixed 
all free parameters, such that the theory is now fully predictive for all other processes that can 
be described by our one-loop HEL. 

QED is usually defined at some high ultraviolet (UV) scale in terms of a bare action. It 
is predictive for electrodynamic processes at low energy scales to a high accuracy, once all 
renormalization-group "relevant" parameters are specified. QED has two such parameters, the 
coupling e and the electron mass m. During the derivation of the HEL (1), we have to give a 
prescription how to fix these parameters to their physical values, since they turn out to be scale- 
dependent. However, since the HEL of Eq. (1) describes processes with external photon lines 
to one loop only, the scale dependence of m remains invisible in this calculation. Hence, to this 
order of accuracy, we can fix the electron mass to its value measurable at the scale of, say, atomic 
physics, and insert this into Eq. (1). In the HEL language, the scale dependence of the mass be- 
comes visible at two-loop order, see, e.g., 30? 32llH^ . This subtlety emerges in the following way: 
The unrenormalized two-loop expression contains a term which may be written as 6m'^ {d/dm?) 
(one- loop Lagrangian), where m is the electron mass and 5m the radiative modification (the 
physical mass is then m^Yi — '^'^ + ^m^ . The original Maxwell Lagrangian is independent of the 
electron mass. The one-loop Lagrangian that we work with depends on the electron mass. The 
two-loop Lagrangian has effective self-energy corrections to the dressed fermion propagators that 
interact with themselves through a further radiative photon. This self-energy effect is known to 
produce a radiative renormalization 5m of the electron mass. Therefore, the above term has to 
be interpreted as a correction to the electron mass that enters into the one-loop Lagrangian, 
expressing the fact that "one-loop Lagrangian +5m'^ {d/dm^) (one-loop Lagrangian) = one-loop 
Lagrangian as a function of "iRenormaiized ~ + (5m^." This is the mechanism by which the 
anomalous mass dimension creeps into the renormalization group analysis of the action [HO] . 

Renormalization— Group Invariance: The QED Effective action is well known to be a mod- 
ification of the Maxwell Lagrangian relevant to the situation of strong background fields. By 
contrast, vacuum polarization modifies the Coulomb at small distances (large momenta). How- 
ever, there is a connection between the two regimes, asymptotically. Specifically, Ritus (see |32j 
and references therein) has shown that an interesting connection exists between the corrections 
to the Maxwell equations in intense fields and vacuum-polarization corrections to QED at large 
momenta. The connection can now be found upon considering Eqs. (25) and (26) of Ref. 32 . Ob- 
serve that the mass renormalization in Eq. (14) of Ritus also has the right asymptotic behaviour 
at small proper time sq . One of the basic ideas of the renormalization group invariance is that the 
product of field and charge, which enters into the covariant coupling, — e^A^^ — > — eA^ in 
going from bare to renormalized quantities, necessarily has to be preserved under the renormal- 
ization because space itself {d^ = d/dx^) does not stretch under the renormalization of charge 

and field. The charge and the field are renormalized as cq = ^/^e and A^ = Z^^^A^. [H ere, 
"0" denotes the bare quantity.] 

Let us finally discuss the scale dependence of the coupling e as it occurs in the HEL. Starting the 
computation with a bare (Maxwell) action at UV scale A with a bare coupling cq = eo(A), we 
obtain a A- and eo-dependent contribution to the Maxwell action ~ J^. Together with the bare 
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action, we can trade these contributions for a scale-dependent coupling e = e(/u) that describes 
the interaction of a photon with an electron at an energy scale /i [31301132] • Of course, /i is an 
arbitrary parameter, and the requirement that no physical process should depend on n defines 
how the value of the coupling parameter e(/x) has to be changed upon a variation of For 
quantitative predictivity, we have to fix the value of e(^) at one particular scale ^. Starting 
from this consideration, one may derive functional relationships for the asymptotic behaviour 
of the QED Effective Action in intense fields, which in turn enable the direct transition to 
renormalization-group equations of the Callan-Symanzik type [see Eq. (45) of ^2]- Other types 
of renormalization-group equations, for example those of the Gell-Mann Low type, including a 
thorough discussion, can be found in 5 and j30j . 

Of course, the dependence on fj, of the HEL could be formally kept such as in JOl which would 
correspond to not specifying a renormalization condition. However, once the input of a single 
experiment (e.g., Thompson scattering) is taken into account, no freedom is left in the param- 
eters, the ^-dependence has disappeared and QED is fully predictive. In particular, there is no 
room for further "logarithmic correction terms" |lUj . 
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